Quantum electric dipoles in Dimer Mott insulator coupled to spin degrees of freedom 
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We present an effective dipolar-spin model based on the strong coupling analysis, which may explain the 
possible origin of "spin liquid insulator". The issue is related to a dimer Mott insulator reminiscent of an 
organic triangular lattice system, k-ET2Cu2(CN)3, whose gapless spin liquid state has been discussed in the 
context of geometrical frustration of exchange coupling, J, between spins on dimer orbitals. It turns out that 
another degrees of freedom within the insulator, quantum electric-dipoles on dimers, interact with each other 
and significantly suppress J through the dipolar-spin coupling, resulting in a possible "dipolar-spin liquid" . 

PACS numbers: 75.40.-s, 72.80.Ng, 72.80.Le, 75.5().Mm, 77.22.Ej 



A well known terminology "spin liquid" continues to at- 
tract attention ever since the Anderson's resonanting valence 
bond picture was presentedfT]. This state of matter is recently 
discussed in organic solid, K-(BEDT-TTF)2Cu2(CN)3||^j5Il 
(BEDT-TTF abbreviated as ET), in the context of absence of 
magnetic orders and of the possibility of gapless excitations. 
The solid has a triangular lattice structure at half-filling in unit 
of dimerized molecules, and is driven towards a Mott insulat- 
ing state by the strong on-dimer Coulomb interaction l'6'l. The 
origin of its spin liquid state is then ascribed to the frustration 
effect from the nearly regular triangular geometry of the spin 
exchange coupling, J^g/ Jcg ~ 1, which is estimated from 
the corresponding effective transfer integrals between dimer 
orbitals, i^g/tcff ~ I.OsItII- Theoretical works are devel- 
oped successively within the half-filled single band Hubbard 
model (SBH) on the triangular lattice|8- 12]; an exotic non- 
magnetic (gapless) insulator is found in the anisotropy range 
of t^g/tcff ^ O.T-lfsl. Some interpretations are given by 
many-body spin exchanges beyond the Heisenberg model and 
the spinon Fermi surface Uslll^jl or by vison excitation llisll . 

Recently, however, k-ET2Cu2(CN)3 is reported to have 
anomaly at T~6K in lattice expansion coefficient! 16] as well 
as in dielectric constant which shows a relaxor ferroelectric- 
like behavior above this temperature fv^. They suggest that 
the electronic degrees of freedom is still active in the insu- 
lating state. Coincidentally, model parameters are replaced 
from the above mentioned ones||3] by the ab-initio calculation 

to W^cff ~ O.SlUS, 19], F/C/dimor -0.4 with t/dimcr/i ^ 

15 111811 . where V and t/dimcr are the inter-dimer and on-dimer 
Coulomb interactions, respectively. Therefore the geometrical 
frustration effect is not strong. Instead, the inter-site Coulomb 
interaction, Vij, is large which shall play certain role in the 
low energy physics of this intriguing state. Since the system 
is quarter-filled in unit of molecule, there is an instability to- 
wards charge order by Vij, which may compete with the Mott 
insulator[20]. In this Letter, we explicitly include Vij for the 
first time in describing k-ET2Cu2(CN)3 and demonstrate an- 
other scenario for the suppression of magnetic orders. We 
describe charge degrees of freedom in the insulating state as 
"quantum electric-dipoles " which fluctuate within the dimer 
Dipoles couple to spins through the inter-dimer fluctuation of 
charges. The suppression of dipolar fluctuation by Vij leads 



to significant decrease of J's and of magnetic correlations. 

We go back to the basic model of organic solids in unit 
of molecule [20], a quarter-filled two-band extended Hubbard 
model, whose Hamiltonian is given as. 
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Here, Cj^ I c-^ are creation/annihilation operators of electrons 



with spin CT(=t, i) and nj^r = 



and rij — Uif + Ui^ are 



number operators. The model includes on-site (U) and nearest 
neighbor (Vij) interactions. We consider strong dimerization 
effect, namely, each pair of sites connected by strong interac- 
tions, {tij, Vij) = {td, Vd), called "dimers" has one electron 
(i.e., half-filled in unit of dimer) on an average. 

The half-filled SBH in Refs.||8-12] is a limiting case of 



quarter-filled Eq.([T]l, td ^ tij and Vij = 0; In Eq.([T]l each 
dimer has sixteen bases which is reduced to four bases in 
SBH by the "dimer approximation". This four bases could ac- 
count for charge fluctuation between dimers and describe both 
metal and dimer Mott insulator. Whereas, details of the two 
dimerized sites are neglected, e.g., a doubly occupied basis of 
SBH, ('|~|) , represents only one bonding state among six dou- 
bly occupied states of Eq.©, (ti, 0), (0, tl), (t, i), (i, t), (t 
, t), a, i), where a/Q denote the presence/absence of spin-cr 
electrons on (site-l,site-2). The intra-dimer charge dispro- 
portionation/fluctuation is not considered, which is no longer 
legitimated when Vij is large (> id)! 18]. 

Instead, we take the strong coupling limit (insulator) of 
Eq.dTJ as tij /U,tij /Vd,Vij /U,Vij /Vd 0. This approxi- 
mation projects out bases having double occupancy on dimers 
with energy of order U or Vd, and describes both dimer 
Mott and charge ordered insulators. The remaining bases 
keep exactly one fermion per dimer as, (site-1, site-2)=(t 
,0), (0,t), (i, 0), and (0,^,]. We introduce electric-dipolar 
and spin operators, P and S. Then, the above four states are 
described in the 2®2 spinor representation as, {P^-,S^') — 
|), (— |, |), (|, -^), and (— The quantization 
z-axis of dipoles is fixed to the dimer-bond direction in real 
space. In the following, we derive the effective Hamiltonian 
by treating tij and Vij (^ Vd) perturbatively, by which the 
charge fluctuation between dimers is taken into account. 
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FIG. 1: (Color online) (a) K-type lattice structure adopted to Eq.(l) in unit of molecule (circle); Bold line represents dimer bonds, (td, Vd)- (b) 
Anisotropic triangular lattice in unit of dimer, transformed from panel (a), (c) Two different types of connections between dimers in panel (a). 
Polarized (d)(e) and unpolarized (f) configurations of quantum dipoles. Charges avoid neighboring alignment of strong Vij shown explictly in 
the fermionic representation (left panels). Panels (g)-(j) are the representative second order perturbation processes, which generate effective 
interactions in Eq.JSj, where filled circle and arrows represent the electrons and spins, respectively. 



For the realization of Eq.([T]l in the bulk system, we 
choose a model lattice of K-type organic solids 12311 shown 
in Fig. [n^); the inter-dimer interactions are {tij,Vij) = 
(tp, V^), {tq, Vq), {ts, Vb) along three different bonds. If we 
take dimer as a unit, this lattice is mapped to an anisotropic 
triangular lattice in Fig. \llh), whose horizontal (W) and di- 
agonal (W) inter-dimer bonds originate from (isiVe) and 
{tp,Vp), {tq,Vq) of Fig. [TJa), respectively. The effective 
Hamiltonian, Ucs = H^^^ +n^'^'^ +'H^^^ +n^'^\ is generated 
in unit of dimer in Fig.[Tlb). The first order Hamiltonian is. 
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where W(5 — Vb and Wq ~ Vq — Vp. Eq.(|2]l is nothing 
but a transverse Ising model of dipoles, where a competi- 
tion of correlation and local quantum fluctuation of dipoles is 
imbedded, namely the Ising interaction term Wq{^ Vij) ver- 
sus the transverse field td- Representative dipolar states are 
shown in Figs. [nd)-[Tlf); when td -C Vij, there are two dif- 
ferent spacial orders depending on the geometry of Vij. At 
td ^ Vij, dipoles fluctuate and stay spacially uniform. The 
former "dipolar-solid" corresponds to charge order and the lat- 
ter "dipolar-liquid" to dimer Mott insulator. At this order, the 
spin degrees of freedom are fully degenerate. 

Degeneracy of spins are lifted at the second order level. The 
second order terms for two bonds (W^' = W' and W 12411 ) 
yields. 
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-wl^ p^p^ + wLSa, ■ - w^^Jp^p^) {s^ ■ s, 



+WliP+ + P-)- Wl^siP^ + p-){Sa. ■ Sp) . (3) 
Let us first focus on bond-T4^'. The first (P^P|) -term origi- 



nates from the process in Fig. [Hg)- Noteworthy is the emer- 
gence of a dipolar-spin coupling term, {P^Pp){Sa ■ Sp), 
which together with VF'ss-term originate from the process 
in Fig. [Ilh); exchange of spins occurs only when fermions 
occupy nearest neighbor sites, i.e., when dipoles are an- 
tifeiToelectric {PajPp) = 5)- This four-body term 
reminds of the Kugel-Khomskii Hamiltonian discussed in 
manganites['2ll'22l] in the context of orbital-spin coupling. 

The diagonal bond-W^ consists of two interactions {p, q), 
each generating terms in the same manner as bond-M^'. In ad- 
dition, there are processes shown in Figs.llJi) and[Tlj); Dimers 
exchange their fermions through two connections p and q by 
ending up flipping P^, namely, fermion in dimer-a changes 
its site location. This generates the last two terms in Eq.(|3]i 
only for bond- VF, i.e., = W^,, = 0. 

Spins follow dipoles through these second-order dipolar- 
spin coupling terms in Eq.(|3]i. Effective interactions between 
spins on neighboring dimers {a, (3) are evaluated from the ex- 
pectation values of coefficient of Sa ■ Sp of Eq.Q as, 

J = - W^pps.(P„^P|) " (M^pss(P+ + P-)). (4) 

In order to examine the actual competition between td and 
Vij , we perform exact diagonalization on "HefF in finite dimer 
clusters 125,1. We adopt {tp, tq) = (0.7, -0.25) in unit of is = 
1 which describes K-ET2Cu2(NCS)3|l7jt], and take {V,Vd) = 
(15, 10). These values are interpreted to coefficients and 
of Hcff. For the choice of Vij, we take > and Vp ~ 
Vb ~ 0, concentrating on the type of dipolar solid given in 
Fig lfe), which is sufficient to clarify the essential physics of 
the competition of Vij and td- 
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FIG. 2: (Color online) Results on Heft as functions of td and Vq 
with fixed {tB,tp,tg) = (1,0.7,-0.25) and ([/, V^s) = 

(15,10,0,0). (a) td-Vq phase diagram. Along its Vq = 3 line, 
(b) J and J', (c) (P<^P|) and (P+) are shown. The broken lines, 
Jd and J^, are the td ^ oo limit of J and J' obtained by taking 
(PqP|) = and (P±) = -0.5 in Eq.©. Solid line, J^fi, is the 
evaluation from the conventional dimer approximation. Panels (d)- 
(f) give the comparison up to second, third, and fourth order pertur- 
bation; (d) J'/ J and the structural factors of (e) spins XQ and (f) 
charges Cq with Q — {n,n). The bulk ground state of Heisenberg 
model of the correpsonding J'/ J is shown on the right of panel (d). 



We first elucidate the phase diagram on the plane of td and 
Vq as shown in Fig. EJa). As anticipated, the dipolar-solid 
(Fig. [Ite)) and liquid (Fig. [TJf)) appears at large Vq and td, 
respectively. The solid-liquid phase boundary is determined 
as a minimum of charge gap of the U — c>o-limit of Eq.([T]l. 

Next, we vary td along the fixed line of = 3 in the phase 
diagram at the second order level, T-LeS = 'H^^-' + H'^-*. Fig- 
ure |2jb) shows that J' is suppresssed in the dipolar-solid at 
small td, and increases significantly by td- This can be ex- 
plained by the remarkable t^-dependence of {P^P^) shown 
in Fig. |2jc); In Eq.(|4|i, the second term of J' with the con- 
stant coefficient, W^^^^ = ^*_y > 0, has large negative con- 
tribution to J' by {P^Pp) =^ 0.25 at ^ 1, which goes to 
zero as {P^Pf) — >^ at — > oo. Similar discussion holds 
for J; the second and third term of Eq.® are both positive 
but are decreasing and increasing functions of td, respectively, 

since Wppss = -^zV 



j^i2_>o,(P„^P|),, -0.25^0, and 

VFpss > 0, {Pa) — — > —0.5. The two variations cancel out, 
keeping J almost unchanged by td- 

One important point is the unexpected decrease of J' due 
to Vji- In general, K , works to screen U and to enhance 



J 112611 . This effect actually appears in the increase of de- 



case Vij works directly on dipoles at the first order level, and 
varies (P^Pp) more significantly than their coefficients, sup- 
pressing J'- By contrast, J^^g = 4t^fj ^/?7dimer obtained by 
the conventional dimer approximation based on SBH shows a 
"screening effect" as shown in Fig.|2jb), where t'^^ — and 

^dimcr — *^td ^ (l 



4td 
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nominator of Wppss oc {U — Vij) ^ hy Vij- However, in our 
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Another interesting aspect is the strong correlation between 
J and J'. In the dipolar-solid state, the dipolar vectors are 
bound to particular spacial directions, which contribute to 
both J and J' through (P^P|) - ±0.25. The antiparallel 
dipoles along bond-iy' always suppresses J', whereas along 
bond-M^, the antiparallel dipoles in Fig. [TJd) suppress J but 
parallel ones in FiglTJe) do not. In both cases, however, we 
confirmed that J'/ J is suppressed at small td (see Fig.|2|d)). 

In this way, at small td, the anisotropy of J'/ J becomes 
large and the geometry of spin interactions approaches a 
square lattice- For reference, we consider the bulk ground 
state of the Heisenberg model] 27 -3^; at J' / J ^ 0.7 the an- 
tiferromagnetic order sets in, which is reflected in the spin 
structural factor, xq, as a single peak at Q = (tt, tt). Fig- 
ure |2je) shows the peak amplitude of XQ of our model ver- 
sus the Heisenberg model calculated on the same finite clus- 
ter The phase transition of the spin sector takes place at 
J' /J - 0.7, i.e., at (trf/t), - 2.4. At td/t > (td/t)^ 
well inside the dipolar-liquid phase, the present model has 
the same XQ the Heisenberg one. By contrast, XQ is 
significantly suppressed from that of the Heisenberg one at 
td/t ^ {td/t)c- This indicates that in approaching the dipolar- 
solid phase, the dimer spin degrees of freedom even at second 
order is no longer mapped onto the Heisenberg model- In such 
case, the multiple spin-exchanges which emerge as higher or- 
der perturbation terms may become important. Therefore, 
we calculate all the perturbation processes up to fourth order 
-Hoff = -^(1) + + _^ -^(4)^ eight-bodWnterac- 
tions, and including the four-dimer ring exchangesyjj] which 
maximally amount to J^/t ~ 0.0027). The numerical results 
up to third and fourth orders are compared with the second or- 
der ones in Fig.|2|d)-(f); J'/ J as well as the structural factor 
of dipoles (Cq in Fig. Eff)) are almost unchanged. By con- 
trast, xq shows a significant upturn {td/t < 1.7) at the fourth 
order level lf32ll . Still, in the vicinity of the dipolar-solid/liquid 
phase boundary, td = 1.7—2.4, there remains a strong suppres- 
sion of antiferromagnetic correlation at all orders. The cor- 
responding region is descibed as "dipolar-spin liquid" phase 
in Fig. 12 a). It apparently not originates from spin degrees of 
freedom alone. We consider that the liquid dipoles develop 
their short-range correlation, which rumple the distribution 
of charges, resulting in a spacially non-uniform value of J's 
within the certain imaginary timescale, leading to the possi- 
ble spin liquid state. Thus, in the dimer system, the four-body 
dipolar-spin interactions in Eq.© play comparably important 
role as a well-known higher order ring exchange terms li31il in 
liquidizing the spin sector. 

To summarize, we derived an effective dipolar-spin model 
relevant in the insulating phases of the quarter-filled extended 
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Hubbard model with lattice dimerization. The intra-dimer 
charge degrees of freedom is described for the first time as 
quantum electric dipoles fluctuating by (the inter-dimer 
transfer integral). At large the conventional dimer Mott 
insulator, namely a dipolar-liquid, is stable. The inter-dimer 
Coulomb interaction, V, competes with td and at large V/td 
a dipolar-solid (charge order) emerge. The spins couple to 
the dipoles through the inter-dimer charge fluctuation, and 
in approaching the dipolar-solid from a liquid phase, strong 
modicication of spin-spin interactions, J and J', as well as 
significant suppression of the antiferromagnetic correlation is 
observed, indicating the existence of a "dipolar-spin liquid". 



We anticipate that the so-called spin liquid Mott insulator 
K-ET2Cu2(NCS)3 lies in our "dipolar-spin liquid" phase at 
low temperature, where both spins and charges remain short 
range ordered. The final comment is given on the com- 
parison of parameters with experiments and other theories. 
K-ET2Cu2(NCS)3 undergoes a metal-insulator (MI) transi- 
tion at ~0.4GPa||2[. The ab-initio calculation shows that t^j 
varies by 1.3 times from ambient pressure to 0.75GPa|19], 
which brings Udima/tcH from 15 to ll lllSi 13311 . This value 
is already larger than the conventional empirical estimate of 
Kiimer/ ^eff ^ 7— 8|j2|] . In such casc, our perturbative treatment 
up to fourth order can be qualitatively adopted, at least at am- 
bient pressure. Starting from the same Uj [7] the conventional 
dimer approximation on a SBH reaches the effective Heisen- 
berg coupling, J^fj / Jeff ^ 1, whereas our dipolar-spin-liquid 
always has J'/ 0.6— 0.7. Interestingly, our J'/ J coincides 
without assumption into a result of the ab-initio calculations 
^ O.siillil]. Thus, we interpret the ab- 



as, y/j^r^t'^ff/Us' 

initio results as a renormalized value of our effective J's after 
including Vij. Regarding the interpretation of 6K anomaly, 
we expect the following senario; at high tempearture (T), the 
thermal flucuation of dipoles is dominant (dipolar liquid is sta- 
ble). At T 0, the dipoles remain short range ordered as well 
due to the quantum fluctuation (td), and so as the spins which 
follow dipoles. In lowering T, the dipolar(electronic) corre- 
lation once grows, but still remains a liquid towards T = 0. 
T — 6K is possibly a maximally correlated point (maximum 
V/ td), where J'/ J is suppressed at most. In fact, a relaxor fer- 
roelectric behavior of dielectric constant indicates a pseudo- 
transition at T^=6K [.17il . The lattice anisotropy also takes a lo- 
cal maximum here, which implies that J' / J also pass through 
an extreme value. To confirm this senario, further develop- 
ment in both theories and experiments are requked. Thus, so 
far, one cannot exclude the possibility that the "gapless spin 
liquid" may not be the result of the geometrical frustration but 
of a strong correlation between spins and charges(dipoles). 
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